Dynamic localization and transport in complex crystals 
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The behavior of a Bloch particle in a complex crystal with VT symmetry subjected to a sinusoidal 
ac force is theoretically investigated. For unbroken VT symmetry and in the single-band approxima- 
tion, it is shown that time reversal symmetry of the ac force preserves the reality of the quasienergy 
spectrum. Like in ordinary crystals, exact band collapse, corresponding to dynamic localization, is 
attained for a sinusoidal band shape. The wave packet dynamics turns out to be deeply modified at 
the VT symmetry breaking point, where band merging occurs and Bragg scattering in the crystal 
becomes highly non-reciprocal. 

PACS numbers: 73. 23. Ad, 11.30.Er, 72.10.Bg 



INTRODUCTION 



The dynamics of matter or classical waves in peri- 
odic potentials subjected to external dc or ac forces 
is strongly influenced by Bragg scattering, which is re- 
sponsible for such important effects as Bloch oscillations 
(BOs) and dynamic localization (DL). Dynamic local- 
ization was originally proposed by Dunlap and Kenkre 
d as a suppression of the broadening of a charged- 
particle wave packet as it moves along a tight-binding 
lattice driven by a sinusoidal ac electric field. DL was 
then explained in terms of quasienergy band collapse [2| , 
and the general conditions of DL beyond the nearest- 
neighboring tight-binding (NNTB) approximation used 
by Dunlap and Kenkre were subsequently investigated 
in Ref. [3] . The interest on DL has been recently renewed 
since the first experimental observations of DL have been 
reported for matter waves trapped in dynamical optical 
lattices |4j6[ , and for light waves in curved waveguide 
arrays |7Hll|. Matter or light waves may also interact 
with complex potentials. Complex crystals for matter 
waves emerge, for instance, in the near resonant inter- 
action of light with open two- level systems [12| , whereas 
in optics complex crystals are realized by waveguide ar- 
rays that include gain and/or loss regions [13]. As com- 
pared to ordinary crystals, complex crystals exhibit some 
unique properties, such as violation of the Friedel's law 
of Bragg scattering and nonreciprocal diffraction [I2I-HH . 
A special class of complex crystals is provided by com- 
plex potentials possessing parity-time (VT) symmetry 
fl3[ [lEl Htjj . An important property of VT crystals is to 
admit of an entirely real- valued energy spectrum below 
a phase transition (symmetry-breaking) point, in spite of 
the non-Hermiticity of the underlyin g H amiltonian [15j . 
A recent study on BOs in VT crystals [17j has shown that 
the common wisdom of coherent quantum transport in a 
crystal is greatly modified when dealing with a complex 
crystal. 

It is the aim of this work to investigate the coherent mo- 
tion of wave packets in a complex VT crystal driven by 
an ac-like force. In particular, it is shown that in the un- 
broken VT symmetry region time reversal symmetry of 



the ac-like force preserves the reality of the quasienergy 
spectrum, and that a full band collapse, corresponding to 
dynamic localization, occurs within the single-band and 
NNTB approximations like in ordinary crystals. How- 
ever, the transport properties of the lattice are deeply 
modified at the VT symmetry breaking, where Bragg 
scattering in the crystal becomes highly non-reciprocal. 
In the following analysis, we will consider specifi cally 
wave packets in a photonic lattice system [7|, LUU, [l3|, UJ| , 
however the results can be applied to other lattice real- 
izations, such as to matter wave tunneling in dynamic 
complex optical lattices. 



II. DYNAMIC LOCALIZATION IN COMPLEX 
CRYSTALS WITH UNBROKEN VT SYMMETRY 

In optics, the coherent motion of charged quantum par- 
ticles in periodic potentials driven by an ac electric field 
can be mimicked by the propagation of light waves in a 
periodically-curved photonic lattice [7|, [l8| . In the scalar 
and paraxial approximations, light propagation at wave- 
length A in the lattice is described by the Schrodinger- 
type wave equation 



iXd z ifj 



2n, 



■d 2 x ^^V{x)^-Fx^ = {UQ-Fx)^. (1) 



where A = X/(2tt) is the reduced wavelength, n s is the 
substrate refractive index, V(x) = n s — n(x) is the po- 
tential, n(x) = n(x + a) is the refractive index profile 
of the lattice (spatial period a), and F = F(z) is a fic- 
titious refractive index gradient proportional to the lo- 
cal waveguide axis curvature which mimics the action of 
a driving force [18|- In particular, a sinusoidal ac-like 
force is simply mimicked by a sinusoidal bending profile 
of the waveguides @- In a complex lattice, the refrac- 
tive index is complex, and the VT symmetry require- 
ment V{— x) = V*(x) corresponds to suitable combina- 
tions of optical gain and loss regions in the lattice as 
discussed in [13]. The real and imaginary parts of the 
potential are denoted by Vr(x) and aVi(x), respectively, 
where a > is a dimensionless parameter that measures 
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the anti-Hermitian strength of Ho- The spectrum of Ho 
turns out to be real- valued for a < a c , where a c > 
corresponds to the transition from unbroken to broken 
VT symmetry. Numerical studies generally show that 
for a < a c the spectrum is composed by bands separated 
by gaps like in an ordinary crystal, whereas for a > a c 
band merging is observed with the appearance of pairs 
of complex-conjugate eigenvalues |13l. Il6j. For instance, 
for the potential defined by 

Vr(x) = V cos(2nx/a) , Vi(x) = V sin(27nr/a), (2) 

one has a c = 1 pjjj]. 

In this section we consider the unbroken symmetry phase, 
i.e. the case a < a c . As recently shown in Ref . , 
for a < a c the motion of a Bloch particle in pres- 
ence of an external dc force F can be described fol- 
lowing the same lines as in ordinary crystals by ex- 
panding the field ip(x, z) as a superposition of Bloch- 
Floquet eigenfunctions <fi n (x,K) = u n (x, k) exp(i^) of 
Ho, where the wave number n varies in the first Bril- 
louin zone, i.e. — fc#/2 < n < fc#/2, ks = 27r/a is the 
Bragg wave number, n is the band index, and u n (x,K) 
is the periodic part of the Bloch-Floquet eigenfunc- 
tion. After setting ip(x,z) — J <iftc n (ft, z)0 n («, z) 
and assuming normalized eigenfunctions such that 

/ dx<j>n>{—X, -K, f )(f) n (x, K) = V n 5 n , n '5(K - K f ) with V n = 

±1, the evolution equations for the spectral coefficients 
c n (ft, z) read [17[ 

iX (f* + ° n = En ^° n ~ FV n^X U) i(K)ci (3) 

where E n {n) is the energy of (j) n (x,K) 
[with E n (-K.) E n {n)\ and X ni i{n) 

(2m /o) f Q dxUn(—x,—K)d K ui(x,K). The off-diagonal 
elements X n? / (n ^ I) in Eq.(3) are responsible for inter- 
band transitions, i.e. Zener tunneling (ZT). If bands n 
and / are separated by a large gap and the ac force F(z) 
is small enough such that \FX u j(k)\ <C \E n (n) — Ei(k)\ 
in the entire Brillouin zone, ZT is negligible as in 
ordinary lattices and one can make the single-band 
approximation by setting X Ui i ~ for n ^ I in Eq.(3). 
In the single-band approximation one thus obtains 

iX (d z + jd K ^j c(z, k) = [E(k) - iF(z)<$>(K)} c(«, z) 

where we omitted, for the sake of simplicity, the band in- 
dex n and set i<&(n) = V n X n ^ n (n). Because of the sym- 
metry of V(x), Re(u n (k,x)) and Im(u n (k,x)) have well 
defined and opposite parity under the inversion x —> — x; 
this implies that is a real- valued function of ft, van- 
ishing for a real potential (i.e. for a = 0). As previ- 
ously shown in [l7j, when a dc force F is applied to the 
crystal, from Eq.(4) it follows that the energy spectrum 
is described by a complex- valued Wannier- Stark ladder. 
The non-reality of the energy spectrum comes from the 
extra-term in Eq.(4) and is physically due to the 



fact that the external dc force F breaks the VT sym- 
metry of the full Hamiltonian H = Ho — Fx. For an 
ac-like force with period A = 27t/cj, because of the z- 
periodicity of the Hamiltonian its energy spectrum is re- 
placed by a quasi-energy spectrum. Moreover, in the 
single-band approximation DL corresponds to a complete 
collapse of the quasi-energy band like in an ordinary crys- 
tal |2[. According to Floquet's theorem of periodic sys- 
tems, the quasi-energy £(k) for the n-th lattice band 
can be readily calculated by looking for a solution to 
Eq.(4) of the form c(z,k) = a(z, k) exp[— i£(n)z/X\ with 
a(z + A, k) = a(z, ft). One obtains 

£{R) = \f W(k')-^(*)*(«')] (5) 

where we have set ft' = ft — fc(A) + k(z) and k(z) 
( 1 / A ) f* d£F(£). Let us assume that the ac forcing F(z) 
is an odd function with respect to some point zo, i.e. 
that F(z — zo) = —F(zo — z) for some zo in the oscil- 
lation cycle. This condition is satisfied, for instance, for 
the important case of an harmonic (e.g. sinusoidal or 
cosinusoidal) ac driving force, originally considered by 
Dunlap and Kenkre [l| and that will be assumed in the 
following. Owing to this additional temporal symmetry 
on the driving force, which is absent for the BO prob- 
lem [17], a real- valued quasi-energy spectrum for the 
non-Hermitian time-periodic Hamiltonian Ho — F(z)x 
is obtained. In fact, in this case the imaginary term 
on the right hand side of Eq.(5) vanishes after integra- 
tion because F(z) and have opposite parity for 
the inversion (z — zo) — » — (z — zq). The quasi-energy 
spectrum is thus real-valued and its expression takes 
the usual form as in a conventional crystal. DL corre- 
sponds to a collapse of the quasienergy band £(ft), i.e. to 
d£(K)/dhz = 0. For most driving fields like a sinusoidal 
field, DL can be attained exactly solely in the NNTB 
approximation, i.e. when the band shape E(k) is sinu- 
soidal [3], E(k) = Eo — Acos(fta). In this case, assum- 
ing for the sake of definiteness a sinusoidal ac-like force 
F(z) = Focos(ujz), the explicit form of the quasi-energy 
reads [2[ 

£(k) =E - AJo (J^j cos(fta). (6) 

Band collapse, leading to DL, is thus attained when 
Jo{Foa/Xuj) =0 0,0,0- 

To check the correctness of the analysis, we investigated 
DL for the complex crystal V(x) = Vb[cos(27nr/a) + 
iasm(27Tx/a)] in the unbroken VT symmetry phase (a < 
1) by a direct numerical analysis of Eq.(l) using a pseudo- 
spectral split-step method with absorbing boundary con- 
ditions. Figure 1 shows a typical example of band dia- 
gram below the phase transition point (a = 0.3), numeri- 
cally computed by a spectral analysis of the Hamiltonian 
Ho- Note that the lowest band of the array is with ex- 
cellent accuracy approximated by a sinusoidal curve and 
turns out to be separated by a large gap from the second 
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FIG. 1: (color online) Band diagram of the complex potential 
denned by Eq.(2) for parameter values A = 633 nm, a = 
8 /m, n s = 1.42, V = 0.002, and a = 0.3. The dotted 
curve shows, for comparison, the parabolic dispersion relation 
E(k) = (Aft) 2 /(2n s ), folded inside the first Brillouin zone, 
corresponding to the critical case a = a c = 1. The picture on 
the right hand side is an enlargement of the lowest band (solid 
curve), fitted by a sinusoidal curve (dotted curve), almost 
overlapped with the solid one. 
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FIG. 2: (color online) Wave packet broadening (discrete 
diffraction) and dynamic localization in the lattice of Fig.l 
for broad beam excitation at normal incidence (left panels) 
and for single site excitation at normal incidence (right pan- 
els). In (a) and (b), the external force is absent (Fo = 0). 
In (c), (d) a sinusoidal force with period A = 1 cm and am- 
plitude Fo = 13.32 /im" 1 , corresponding to T = 1.684, is 
applied. In (e) and (f), the forcing amplitude is increased 
to Fo = 19.03 /im -1 , corresponding to the DL condition 
T 2.405. 



band. DL is thus expected to occur provided that the 
lattice is excited in its lowest order band. Figures 2(a) 
and (b) show a typical spreading (discrete diffraction) of 
a Gaussian wave packet ^(x, 0) = exp(— x 2 /w 2 ) in the 
absence of the external force, for either a broad Gaus- 
sian beam that excites at normal incidence a few lattice 



sites at z = [2w = 5a, Fig. 2(a)], and a narrow Gaussian 
beam that excites at normal incidence a single well of the 
lattice [2w = 2a/3, Fig. 2(b)]. In both cases, the lowest 
band of the array is mainly excited, as discussed e.g. in 
After application of the sinusoidal force [Figs.2(c-f)], 
suppression of beam diffraction and self-imaging effects 
are clearly observed when period and amplitude of forc- 
ing satisfy the condition T = Foa/(Xuj) = 2.405 (first 
zero of Bessel Jq function), as shown in Figs. 2(e) and (f). 



III. WAVE PACKET DYNAMICS AT THE VT 
SYMMETRY-BREAKING POINT 

As a is increased to reach and cross the VT symmetry- 
breaking point a c , gap narrowing till band merging, as- 
sociated to the appearance of pairs of complex-conjugate 
eigenvalues, is observed [l3j |. For instance, for the po- 
tential defined by Eq.(2), at the transition point the 
band diagram is given by the free-particle energy dis- 
persion curve E — X 2 k 2 /(2n s ), periodically folded in- 
side the first Brillouin zone [13j (see the dotted curve 
in Fig.l). In this case, as previously noticed for BOs 
[l7| , wave packet transport is deeply modified and can 
not be described by means of the canonical model (3) in- 
troduced in the previous section. From a physical view- 
point, this is related to the highly non-reciprocal behav- 
ior of Bragg scattering in the crystal and violation of 
Friedel's law of Bragg diffraction for crystal inversion 
[l2l [lij . As in Ref. [I?}, we limit here to consider 
the dynamical behavior of a broad wave packet in the 
complex lattice denned by Eq.(2) at a = a c = 1, i.e. 
V(x) = Vo exp(z/c^x), which enables a rather simple ana- 
lytical and physical analysis. As opposed to Ref. [17j, we 
assume here an ac-like force, namely F(z) = Focos(uz). 
Figure 3 shows a few typical examples of wave packet 
dynamics as obtained by a numerical analysis of Eq.(l) 
when the lattice is excited at z = with a broad Gaus- 
sian beam ^(x,0) = exp(— x 2 /w 2 ) for a fixed value of 
the ac modulation period A = 2tt/uj and for increas- 
ing values of the amplitude Fo. For a small amplitude 
Fo, it turns out that the wave packet propagates as if 
the lattice were absent [Fig. 3(a)], following an oscillatory 
path as predicted by the semiclassical analysis of Eq.(l) 
with V = 0. Owing to the external force, the mean 
wave packet momentum varies periodically according to 
Xk(z) = Jq d^F(^) = (Fq/uj) sin(cjz) . As the forcing 
Fo reaches and crosses the critical value F c = /cbAcj/2, 
new wave packets, arising from first-order Bragg diffrac- 
tion, periodically bifurcate from the primary beam at 
propagation distances z satisfying the Bragg condition 
k(z) = —ks/2 [see Figs. 3(b) and 3(c)]. The mean mo- 
mentum of the first-order diffracted wave packets differs 
from that of the primary wave packet by an additional 
term which explains the refraction of the first-order 
diffracted beams at the angle 6 = dx/dz = kBX/n s ob- 
served in Figs. 3(b) and 3(c). At stronger forcing, namely 
for F > 3F C , additional wave packets bifurcate from the 
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first-order wave packets at propagation distances z such 
that k(z) = — 3fc#/2, as shown in Fig. 3(d). These wave 
packets originate from second-order Bragg diffraction in 
the crystal and greatly complicate the pattern scenario. 
The mean momentum of second-order diffracted wave 
packets differs from that of the primary wave packet 
by the additional term 2k which explains the larger 
(twice) refraction angle of second-order diffracted beams 
as compared to that of first-order diffracted beams [see 
Fig. 3(d)]. At even higher forcing, i.e. at F > 5F C , 
new wave packets bifurcate from the second-order wave 
packets because of third-order Bragg diffraction at prop- 
agation distances z such that k(z) = — 5/c#/2, and so 
on. The appearance of wave packets generated by Bragg 
diffraction at various orders is rather abrupt, as indi- 
cated by the behavior of the normalized beam power 
P(z) = J dx\ip(x,z)\ 2 / J dx\i/j(x,0)\ 2 versus propagation 
distance z shown in the right panels of Fig. 3. 
The dynamical scenario observed in numerical simula- 
tions can be analytically captured by considering the lim- 
iting case of a plane wave exciting the crystal at z = 
with initial wave number k = 0. In fact, the solution to 
Eq.(l) with the initial condition ip(x,0) = 1 is given by 
the superposition of diffracted plane waves at different 
Bragg orders according to 

oo 

i/j(x, z) = ^ a n (z) ex.p[ik(z)x + inksx — ij n (z)]. (7) 



In Eq.(7) we have set 



k(z) 



7n0) 



Xoj 



sm(uoz) 



^- [ Z d^[nko + k(0} 2 , 
2n s J Q 



(8) 
(9) 



whereas the amplitudes a n (z) are calculated from the 
recurrence relations 



a n (z) = -iy J dfa n _i(f)exp [ij n (z) - Hn 



(10) 

with ao(z) = 1. The amplitude ao corresponds, in Fig. 3, 
to the primary wave packet, and the independence of ao 
from z indicates that this wave packet propagates as if 
the lattice were absent. The amplitude a\ corresponds to 
the first-order diffracted wave packets, ci2 to second-order 
diffracted wave packets, and so on. According to Eq.(10), 
a wave packet corresponding to Bragg diffraction at or- 
der n bifurcates from a wave packet of order (n — 1), 
and Bragg diffraction is effective provided that the phase 
difference ip n (z) = j n (z) ~ ln-i( z ) entering in the expo- 
nential of the integral on the right-hand-side of Eq.(10) 
has a stationary point. The condition d(p n (z)/dz = is 
satisfied at propagation distances zq such that 



k(z ) 



-k f 



I) (n = 1,2,3,...). (11) 




FIG. 3: (color online) Propagation of a broad Gaussian wave 
packet (input beam size w = 80 fim) in the complex lat- 
tice V(x) = Vb exp(27rix/a) for Vb = 0.0002, a = 6 /xm, 
A = 633 nm, n s = 1.42, subjected to a sinusoidal ac-like 
force F(z) — Fo cos(27T2;/A) with period A = 1 cm and with 
increasing values of force amplitude Fo: (a) Fo/F c = 0.5, (b) 
Fo/F c = 1, (c) Fo/F c = 1.5, and (d) F /F c = 3, where the 
critical forcing F c is defined by Eq.(15). Left panels show 
snapshots of \ip(x,z)\; the central panels depict the behavior 
of k(z), normalized to the Bragg wave number fcs, as given 
by Eq.(8); the right panels show the evolution of normalized 
beam power P(z). The horizontal dotted lines indicate the 
crossing points z = zo where new wave packets are generated 
by Bragg scattering. For the sake of clarity, in (d) only the 
first two crossing points are indicated. At the crossing points 
a rather abrupt increase of beam power P(z) is observed. The 
crossing is linear in (b) and (d), and parabolic in (c). 



i.e. for a wave number that reaches the edge of the Bril- 
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louin zone. Note, however, that in Eqs.(7) and (11) n 
is a positive (but not a negative) integer number: this 
circumstance is a clear signature that Bragg scattering 
is highly non-reciprocal, a feature which is peculiar to 
the complex nature of the crystal [HI, [l4[ . If the cross- 
ings of the stationary points are fast enough, a n (z) is 
basically constant far from the stationary points, with 
abrupt changes at z = zq, namely 

a n(^0~ ) - CLn(Zo) + #a n _i(z ) (12) 

where 

fl=-iy J d£exp[zy> n (0] (13) 

and the integral on the right hand side of Eq.(13) is ex- 
tended to the neighborhood of Zq. From these results 
one can readily explain the abrupt changes of the total 
beam power P(z) observed in Fig. 3 (right panels) at the 
crossing points shown in the central panels of Fig. 3, the 
existence of a critical forcing F c below which no bifur- 
cating wave packets appear, as well as the occurrence of 
higher-order bifurcating wave packets and a dynamical 
scenario with increasing complexity as the amplitude of 
forcing is increased. In fact, according to Eqs. (8) and 
(11) the stationary points which generate the bifurcating 
wave packets of order n are obtained from the equation 

£sin(^ ) = -^(2n-l), (14) 

which can be satisfied provided that the amplitude Fq of 
forcing is larger than F c (2n — 1), where 

F c = ^Xujk B (15) 

is the critical forcing amplitude. Therefore, for Fq < F c 
no diffracted wave packets of any order are generated 
[Fig. 3(a)], for F c < Fq < 3F C first-order diffracted wave 
packets are generated, for 3F C < Fq < 5F C first-order 
and second-order wave packets are generated, and so on. 
The graphical determination of the stationary points is 
depicted in the central panels of Fig. 3 as the crossing 
of the sinusoidal curve k(z)/kB with the vertical dashed 
lines k/ks = —0.5 (for first-order Bragg diffraction) and 
k/ks = —1.5 [for second-order Bragg diffraction, shown 
solely in Fig. 3(d)]. Let ut discuss in some detail the gen- 
eration of first-order diffracted wave packets (n = 1). 
For F c < Fo < 3F C [Fig. 3(c)], in each ac oscillation cycle 
there are two stationary points, and the crossing is linear. 
Near each of the stationary points, ipi(z) can be thus ap- 
proximated as (fi(z) = (pi(zo) + (l/2)(d 2 (pi/dz 2 )(z — zo) 2 , 
where the derivative (d 2 (pi/dz 2 ) is calculated at the sta- 
tionary point z = zq. From Eq.(13), it follows that the 
amplitude factor \R\ of the generated diffracted beam 
after each stationary phase point can be approximately 



computed as 




In deriving Eq.(16), we have taken into account that 
(dVi/dz 2 ) = F(z )k B /n s , \F(z )\ = (F 2 - F*) 1 ' 2 and 
d£ exp(z£ 2 ) = For parameter values used in the 
simulations of Fig. 3(c), from Eq.(16) one has \R\ ~ 0.95, 
which is in excellent agreement with the staircase be- 
havior of beam power P(z) shown in the right panel 
of Fig. 3(c). Equation (16) fails to predict the correct 
amplitude factor R when Fo — » F+, i.e. when the two 
crossing points in each ac oscillation cycle coalesce. For 
Fo = F c , the crossing is parabolic [see Fig. 3(b), central 
panel], and <fi(z) in Eq.(13) should now be approximated 
as ipi(z) = (pi(zo) + (l/6)(d 3 (pi/dz 3 )(z — zq) 3 , where the 
derivative (d 3 (pi/dz 3 ) is calculated at the crossing point 
z = zq. In this case, in place of Eq.(16) one has 




where Ai(£) is the Airy function. In deriving 
Eq.(17), we have taken into account that (d 3 (p±/dz 3 ) = 
{ks /n 3 ){dF/ f dz) = \k\uJ 2 j \2n s ) at z = zq, and 
j^o d£exp(i£ s /3) = 27rAi(0). For parameter values used 
in the simulations of Fig. 3(b), from Eq.(17) one has 
\R\ ~ 2.245. Note that, with this amplitude factor, the 
staircase behavior of beam power P(z) shown in the right 
panel of Fig. 3(b) is reproduced with excellent accuracy. 

IV. CONCLUSIONS 

The coherent motion of a Bloch wave packet in a tight- 
binding lattice driven by an ac electric field is known to 
show a self-imaging effect that arises from quasi-energy 
band collapse of the time-periodic Hamiltonian. This 
phenomenon, referred to as dynamic localization [H, 0], 
has been recently observed for both matter and optical 
waves as a suppression of wave packet broadening in the 
lattice. In this work, we investigated theoretically the be- 
havior of a Bloch particle in a complex crystal with VT 
symmetry subjected to a sinusoidal ac-like force. As com- 
pared to ordinary crystals, complex crystals exhibit some 
unique properties, such as violation of the Friedel's law 
of Bragg scattering and nonreciprocal diffraction. For 
an unbroken VT symmetry and in the single-band ap- 
proximation, it has been shown that the quasi-energy 
spectrum of the time-periodic non-Hermitian Hamilto- 
nian remains real- valued. In this regime, like in an or- 
dinary crystal exact band collapse is possible within the 
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NNTB approximation, i.e. for a sinusoidal band shape. 
At the VT symmetry breaking transition point, band 
merging greatly modifies the wave packet dynamics as 
compared to an ordinary crystal. Here we have inves- 
tigated in details the dynamics of a broad wave packet 
in the T^T-symmetric potential V(x) = Vo exp(2irix/ a) 
and shown that the complexity of the dynamical scenario 
greatly increases as the strength of forcing is increased, 
with the appearance of a cascading of bifurcating wave 



packets. The main features observed in numerical sim- 
ulations have been analytically explained by considering 
the Bragg scattering of a plane wave in the complex lat- 
tice subjected to the sinusoidal ac driving force. These 
results are complementary to the recent study of Bloch 
oscillations in complex crystals [17j, and are expected to 
motivate further investigations aimed to explore the ex- 
otic transport properties o complex crystals. 
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